On the effect of secondary protons on baryon and proton number cumulants in 

event-by-event analysis 
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We investigate the effects of secondary (knockout) protons, which constitute about 20% of the 
observed protons at STAR, on the higher order cumulants of proton and baryon numbers measured 
by event-by-event analyses in relativistic heavy ion collisions. We argue that the contribution of 
this background effect on the cumulants is expressed by a simple formula, and that hence their 
effects can be removed in the experimental analysis, ft is discussed that this background effect has 
non-negligible contribution to recently observed proton number cumulants at STAR, especially the 
third-order one, and that the removal of this effect is crucial to investigate the thermodynamical 
properties of the primordial hot medium appropriately. 



I. INTRODUCTION 

The beam energy scan (BES) program at the Relativis- 
tic Heavy Ion Collider (RHIC) has recently accumulated 
an abundance of experimental data on heavy ion colli- 
sions at various collision energies per nucleon, -^snn 0~ 
4] . Because hot media created by collisions with different 
Jsnn follow different trajectories on the temperature (T) 
and baryon chemical potential (hb) plane during their 
time evolution, careful inspection of these data will en- 
able us to understand the global phase structure of QCD 
on the T-^jlb plane. 

Much attention has been paid to fluctuations of con- 
served charges as experimental observables to probe the 
primordial thermal properties of the hot medium 
Fluctuations are characterized by cumulants. It is theo- 
retically pointed out that the cumulants have character- 
istic behaviors on the T-/j,b plane, reflecting the phase 
structure of QCD \<MlM- 

These behaviors can in prin- 
ciple be measured in heavy ion experiments by observ- 
ing fluctuations in event-by-event analyses. Analyses 
on higher order cumulants are being carried out in the 
BES program [3-H, HH • Because cumulants of conserved 
charges can also be measured on the lattice fluc- 
tuations of conserved charges can be used to compare 
these real and virtual experiments as well as theoretical 
predictions. 

For conserved charges, the time-evolution of fluctua- 
tions in a volume V is slow for large V, and hence fluc- 
tuations in the final state reflect those generated in the 
primordial stage when the rapidity range of detection 
is taken sufficiently large @, H, [HI ■ Indeed, the recent 
experimental result on electric charge fluctuation by AL- 
ICE collaboration at the Large Hadron Collider (LHC) 
suggests that the fluctuation created in the quark-gluon 
phase well survives until the final state at the LHC en- 
ergy ^iiiT = 2.76 TeV 0. 
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The ^snn dependence of net-proton number cumu- 
lants has been analyzed by STAR Collaboration at RHIC 
last two years [l|-|3|- It should be first emphasized that 
the net-proton number is not a conserved charge, and 
should not be considered as a proxy of the net-baryon 
number. Whereas the results on the proton number cu- 
mulants are sometimes identified and compared with the 
baryon number ones in the literature, this identification 
could lead one to a wrong conclusion [TH, To ap- 

propriately compare the experimental results with theo- 
retical predictions on baryon number cumulants, the for- 
mer has to be translated into the latter in the first place 

mm- 

Nevertheless, although with such reservations, the ex- 
perimental results on the net-proton number cumulants 
show some interesting features. To interpret the experi- 
mental results on the cumulants, they are often compared 
with the hadron resonance gas (HRG) model [2(| ■ If the 
fluctuations are generated in a thermal medium described 
well by the free hadronic degrees of freedoms, the exper- 
imental results should well reproduce the HRG model. 
Therefore, the HRG model is used as the baseline from 
which the deviation of the cumulants encodes information 
on the non-hadronic and/or non-thermal feature of the 
primordial hot medium. In the HRG model, the probabil- 
ity distribution of the number of each baryon in a phase 
space is Poissonian (the net number distribution is thus 
given by the Skellam distribution) [201 ] . In the experimen- 
tal results on proton number cumulants at STAR, small 
but statistically-significant deviations from this baseline 
are observed [3J. These deviations are expected to be 
due to non-trivial properties of the hot medium from the 
above argument, and theoretical and experimental verifi- 
cations and investigations of these deviations are impor- 
tant for extracting the global nature of the QCD phase 
structure from these experimental data. 

In experiments, various backgrounds affect measure- 
ments of observables. In the measurement of fluctua- 
tion observables in event-by-event analyses, background 
effects would make the distributions of the fluctuations 
approach Poissionian ones. The observation of the devi- 
ation of the distributions from the Poissonian ones thus 
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would be subject to such background effects. In order to 
see the signals of the non-hadronic property of the hot 
medium encoded as the deviation from the Poisson dis- 
tribution, therefore, careful treatment on the background 
effect is necessary. 

In the present study, we focus on the background ef- 
fect of secondary (knockout) protons, which are protons 
knocked out from materials in the detector by reactions 
with particles created in the hot medium. We derive re- 
lations to remove this background effect from higher or- 
der cumulants, and argue that this background effect has 
a non-negligible contribution to the experimental mea- 
surement of the cumulants of net-proton and net-baryon 
numbers. 



II. SECONDARY PROTONS 

A. General properties 

In the present study, we consider the effects of sec- 
ondary protons on the cumulants of the net-proton and 
net-baryon numbers in relativistic heavy ion collisions. 
Secondary protons are protons which come from not fire- 
balls created by heavy-ion collisions, but materials in the 
detector. Protons in the detector can be knocked out 
by reactions with particles from the fireballs, and enter 
the detector. A major part of the secondary protons can 
be identified and eliminated in the experimental analysis 
by introducing a cut for the distance between the recon- 
structed trajectory of the proton and the collision point, 
which is usually represented by d ca in the literature [2l[ . 
Some of them, however, have reconstructed trajectories 
which pass near the collision point, and are misidentificd 
as true protons created in the fireballs. The abundance 
of these misidentified protons for various transverse mo- 
mentum (pt) bins estimated by an event generator for 
STAR for 62 GeV AuAu collision is shown in Table I in 
Ref. [2l| . With this table and the pr-dependence of the 
proton yield, it is estimated that about 20% of protons 
detected in this experiment are the secondary ones for 
the typical px range used in the event- by- event analy- 
sis of proton number cumulants, 0.4 < px < 0.8GeV. In 
the following, we denote the number of misidentificd sec- 
ondary protons in each event as A^ nd . Notice that the 
anti-proton number is not affected by this background 
effect, because materials in the detector do not contain 
anti-particles. 

The background effect by secondary protons is usu- 
ally subtracted from experimental results on average val- 
ues using estimates with event generators. Their effects 
on the second- and higher-order cumulants, however, are 
nontrivial, and to the best of the authors' knowledge, this 
effect has not been taken care of in the existing experi- 
mental results on fluctuation observables. The purpose of 
the present study is to quantify their effects and make it 
possible to eliminate them in the experimental analyses. 

In order to obtain such a relation, let us clarify sev- 



eral properties of secondary protons. First, we assert 
that the probability distribution of the number of sec- 
ondary protons in each event, P2nd(Np nd ), is Poissonian, 

7 , 2„d(A r p nd ) = A Ar p nd e- A /A2 nd !, to a good approxima- 
tion, with A = (Np lld ). To support this statement, we 
note that the probability r that a particle generated in 
the hot medium creates a misidentified secondary proton 
is small, r« 1. In fact, whereas particles of order 10 3 are 
created in each relativistic heavy ion collision, the value 
of ( A^ nd ) is of order 10 [H[ . This shows that r is of order 
10~ 2 , which is significantly smaller than unity. Moreover, 
events where secondary protons are created are uncorre- 
cted with each other, because the knockout reactions 
take place independently in the bulk detector. These 
two properties are sufficient to guarantee the Possionian 
nature of V2nd{Np nd ). Whereas the total multiplicity 
for each event is not fixed even after fixing the collision 
geometry, and the probability to create a secondary pro- 
ton depends on particle species, the Poissonian nature is 
not altered by these features. The Poisson distribution 
is specified by a single parameter A = (Np nd ), and all 
cumulants take a unique value, 

((SN 2 / d r) c = (N 2 / d ), (1) 

where ((SNp nd ) n ) c is the n-th order cumulant of the 
misidentified secondary proton number. 

Next, let us consider the correlation between Np nd 
and the genuine proton number emitted from the hot 
medium, N™ . Because of the creation mechanism of 
secondary protons, Np nd for each event has a strong pos- 
itive correlation with the total multiplicity 7V tot . Since 
the total multiplicity A tot , which is dominated by pions, 
is correlated with iV™ 1 , in general, Np nd and Np° l are also 
correlated with each other. On the other hand, when col- 
lision events with a fixed N tot are concerned, provided 
that r does not depend on particle species or momen- 
tum, the distribution of Np nd , which is Poissonian, is 
completely determined and hence has no correlation with 
Np° l or any other particle numbers emitted from the hot 
medium. In experiments, the number of charged parti- 
cles, A^h, is used as a proxy of A tot , and analyses are 
performed for collision events with finite ranges of iV c h. 
For example, N c h is used as the centrality selection at 
STAR, and central to peripheral collisions are us ually 
classified centrality bins with the size of 5 or 10% hj- 
In the following, we assume that the cut of N c h in the 
experimental analysis is sufficiently narrow to suppress 
correlations of Np nd with N£° , as well as with the num- 
bers of anti-protons, baryons, and anti-baryons. 

In the analysis of cumulants in STAR experiments, cu- 
mulants have been estimated by the weighted method 
[Tij . In this method, cumulants are first determined with 
narrow ranges of iV c h, and then they are averaged with a 
certain wider range of N^. In this method, the correla- 
tion between Np nd with iV™ 1 vanishes almost completely 
when the value of (Np nd ) determined for each range of 
N c h is used. 
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B. Proton number cumulants 



Neglecting the correlation between N£ and Np , one 
can write the probability distribution of the total proton 
number detected by the detector, A£ x p = N™ 1 + N$ nd , 



x S 



(2) 



Then the moment and cumulant generating functions of 
Np Xp are given by 

G cxp (8) = ^P cxp (iV) e Ar0 = G col (9)G 2nd (9), (3) 

JV 



and 



K exp (0) = log Gexp (9) = if co] (6) + K 2nd {9), (4) 



respectively, where the generating functions for N£ and 



Np nd are defined, respectively, as 



G co i(0) 

G 2nd {9) 



N 
N 



exp^ col (0), (5) 
= expK 2nd (9). (6) 



In Eq. (j4j, K exp (6) is given by the sum of the two inde- 
pendent generating functions K co \{9) and K 2nd (9), owing 
to the independence of the two probabilities V co \{N) and 
V-z^diN). Cumulants ((SN) n ) c are defined by derivatives 
of the cumulant generating function as 



d n K 



86 n 



(7) 



9=0 



Cumulants up to fourth order satisfy (SN) C — (N), 
((SN) 2 ) C = ((SN) 2 }, ((6Nf) c = ((SN) 3 ), and 
((SN) 4 ) C = ({SN) 4 ) - 3((SN) 2 ) 2 . From Eq. Q, it is 
clear that the cumulant of total proton number is given 
by the sum of the two cumulants, 



((5N^) n ) c = ((5K° l ) n )c + ((SN 2nd ) 



= ((SN;° l r) c + (N 2nd ), 



(8) 



where in the last equality we have used Eq. (JTJ) . 

Eq. ([U) relates the experimentally observed proton 
number cumulants, ((SNp Xp ) n ) c , with the genuine pro- 
ton number cumulants in the final state in heavy ion col- 
lisions, ((5N™ l ) n )c- With the value of (N 2nd ) estimated 
in experiments, this relation enables one to eliminate the 
effect of secondary protons from higher-order cumulants. 



C. Net-proton and baryon number cumulants 

Next, we apply a similar argument to the cumulants of 
net-proton and net-baryon numbers, N£_p and A^_g, 



respectively. With the same argument addressed in 
Sec. Ill At one can regard that the correlations between 
N 2nd and the number of each baryon is well suppressed 
with a sufficiently narrow N c h bin. Assuming the absence 
of the correlation in a set of collision events analyzed, the 
probability distribution to observe N^ xp protons and Np 
anti-protons in the experiment is written as 

G c , P (N^,Np) = E Gcoi(N c p °\N p )V 2nd (N 2nd ) 

X <^ XP , ATcol + AT^nd, (9) 

where Gcoi(Np° l , Np) is the probability distribution func- 
tion that the final state of each collision event has A^ 01 
protons and N p anti-protons in the phase space cov- 
ered by the detector. Note that the anti-proton number 
Np does not receive modification from secondary pro- 
tons, and the distinction between "exp" and "col" is not 
needed. This relation immediately leads to 



((SN;^) n ) c = ( (6Np p ) n ) c + (AO- 



(10) 



In Refs. [TaQjl], formulas to obtain net-baryon number 
cumulants using experimentally-observed proton number 
fluctuations were derived. These formulas are based on 
the fact that the probability to observe a baryon coming 
from the hot medium as a proton or a neutron in the 
detector is even and uncorrelated for each baryon to a 
good approximation. This property leads to 

g c oi(N p ,Np)= E •FcoiCNb.JVb) 

x B 1/2 (N p ;N B )B 1/2 (N p ;N & ), (11) 

where J"coi(iVB) Ag) is the probability distribution to 
have A^b baryons and Afg anti-baryons in the phase space 
in the final state, with the binomial distribution function 



B r (N; M) =r JV (l-r) 



\M—N 



Ml 



Nl(M — N)\ 



(12) 



Neglecting the correlation between (anti-)baryon num- 
ber and A" 2nd , one can substitute Eq. ([IT]) into Eq. ((9|), 
obtaining 

(( < 5AT-P g) n ) (rec) = {{S N<%) n ) c + 2(N 2 » d ), (13) 

where ((<5A r g XP g)")i rcc ' 1 is the n-tb. order cumulant of the 
net-baryon number reconstructed from experimentally- 
observed proton number fluctuations using Eqs. (9)-(12) 
in Ref. [18j, while ((SN^_Q) n ) c is the genuine net-baryon 
number cumulants in the final state. 

By employing Eqs. dTUJ) and (fTB")) . one can remove the 
effects of the secondary protons from the cumulants of 
net-proton and net-baryon numbers, respectively. 
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D. Effects on ratios of cumulants 

Next, let us consider how the recent experimental re- 
sults on the net-proton number cumulants are modified 
with the effect of secondary protons using Eq. (ITU|) . 

Because each cumulant is an extensive quantity and 
proportional to the volume of the system, experimental 
results on proton number cumulants are usually discussed 
in terms of their ratios in order to cancel out the volume 
dependence. Eq. (fTU)) shows that the effect of secondary 
protons modifies these ratios as 



which would not be negligible in the discussion of non- 
hadronic nature of the primordial fireballs. 

While we have discussed the effect of secondary pro- 
tons on net-proton number cumulants in this subsection, 
when one wants to compare the experimental results with 
theoretical predictions, it is important to obtain the net- 
baryon number cumulants [l8t Il9| . When the baryon 
number cumulants are analyzed in experiments, the effect 
of secondary protons should be removed with Eq. (fT5|) , 



R 



-COl 

nm 



col 



exp )n> 



2nd\ 



( )™> c ((SN;^) C - (jV-2nd) ■ 



(14) 



This relation shows that when the ratio is unity, i.e. 
the Poissonian value, the effect of secondary protons 
does not alter the ratio. When the ratio = 
({SN^ p p ) n ) c /((5Np > p ) m ) c has a deviation from unity, 
however, removing the effect of secondary protons makes 
the difference large, 



|i-iC™l>|i 



i? cxp l 



(15) 



Substituting (A^ nd ) ~ 0.2(A^ x p) in Eq. flU]), one has 



^ ({SN^Y% - 0.2(iV-P) 



Dcxp 



= R%*-r(l-R%*) + 0(r 2 ), 



(16) 
(17) 



with r = 0.2(A^ x P)/((<5A^ xp ) m ). In the last line, r is 
regarded small and higher order terms in r are neglected; 
in the HRG model, r approaches 0.1 for large ^/snn, 
while r ~ 0.2 for ^/snn ^ 20 GeV for even m. 

At STAR experiment, R^ ~ 0.2 at ySii' ~ 200 GeV 
Q. Substituting this value into Eq. p7|) . one has 



ncol 

/too 



R cxp 
^32 



0.1, 



(18) 



with r ~ 0.1. The third order cumulant at large ^/snn 
thus will receive about 50% modification from the effect 
of secondary protons. 

The values of i?4 Xp at STAR, on the other hand, 
is slightly smaller than but not far from unity for all 
-y/.SNN analyzed in the BES program Q; R^ p ~ 0.8 at 
^/snn — 20 GeV is the smallest value. With this value of 
i?42 P , the difference between i?42 P and R$ expected from 
Eq. (p~7|) is a few percent and is not significant. While the 
correction is not large, we note that the important signal 
is in the difference from the Poissonian value, 1 — R^ . 
In terms of this difference, the effect of secondary proton 
gives rise to about 10% modification for i?42 P — 0.8, 



1 



Rife 1.1(1 



(19) 



III. DISCUSSION AND SUMMARY 



While we have concentrated on the background effects 
by secondary protons on higher-order cumulants of the 
net-proton and net-baryon numbers in this study, similar 
arguments are also applied to other background effects. 
In particular, when a background effect can be regarded 
independent of the baryon numbers in the final state, the 
background effect can be removed in a way completely 
similar to the argument in the present study. One of such 
candidates is the effect of particle misidentification. On 
the other hand, the effects of efficiency and acceptance 
of the detector are a ppr oximately expressed as a bino- 
mial distribution [l8l Il9l. l22l|. and different treatment is 
required [l9j . 

In the present study, we considered the effects of sec- 
ondary protons on higher order cumulants of the net- 
proton and net-baryon numbers measured by event-by- 
event analyses. We argued that the number of secondary 
protons in each event follows Poisson distribution, and 
is uncorrelated with the numbers of each baryon emit- 
ted from fireballs. Using these properties, effects of sec- 
ondary protons on net-proton number cumulants were 
given in a simple form, Eq. (|10p. This formula enables 
one to subtract the background effect from the experi- 
mental results on the cumulants. The formula for the 
net-baryon number cumulants was also given by Eq. ([Tc 



It is found that the effect of secondary protons makes 
the ratios of the cumulants in experimental measure- 
ments closer to unity compared to the ratios with the 
number of the genuine protons coming from the hot 
medium. Such an effect will particularly modify the value 
of the third-order cumulant at large -^/snn- Since the 
deviation of the ratio from unity carries information on 
the non-hadronic and non-thermal properties of the pri- 
mordial fireballs, appropriate treatment on this effect is 
important to reveal the genuine thermodynamical nature 
of the primordial fireballs. 
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